Skyrmion Burst and Multiple Quantum Walk 
in Thin Ferromagnetic Films 
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A giant Skyrmion collapses to a singular point by emitting spin waves in a thin ferromagnetic film, when 
external mangetic field is increased beyond the critical one. The remnant is a single-spin flipped (SSF) point. 
The SSF point has a quantum diffusion dynamics governed by the Heisenberg model. We determine its time 
evolution and show the diffusion process is a continueous-time quantum walk. We also analyze an interference 
of two SSF points after two Skyrmion bursts. Quantum walks for S — 1/2 and 1 are exact solvable. The system 
presents a new type of quantum walk for S > 1 /2, where a SSF point breaks into 2S quantum walkers. It 
is interesting that we can create quantum walkers experimentally at any points in a magnetic thin film, first by 
creating Skyrmions sequentially and then by letting them collapse simultaneously. 



Skyrmions are solitons in a nonlinear field theory, playing 
essential roles in almost all branches of physics[l]. In par- 
ticular, magnetic thin films have recently attracted much at- 
tention owing to real-space observations of Skyrmions y|- A 
Skyrmion crystal^ as well as a single SkyrmionJ^I have been 
identified in chiral magnetic thin films. In spite of their sta- 
bility guaranteed topologically, an intriguing feasibility arises 
that we are able to create them and destroy them experimen- 
tally by breaking the continuity of the field. This is indeed 
the case for giant Skyrmions in magnetic thin films [4]. By 
applying femtosecond optical pulse irradiation focused on a 
micrometer spot, it is possible to destroy the magnetic order 
locally[5]. Then, the dipole-dipole interaction (DDI) gener- 
ates an effective magnetic field, leading to a new magnetic 
order, that is a Skyrmion spin texture||4J. On the other hand, 
as the magnetic field increases beyond the critical one, the 
Skyrmion radius decreases and suddenly shrinks to zero by 
emitting spin waves. This is the Skyrmion burst. 

In this paper we investigate the fate of a Skyrmion after 
its burst. The notion of the topological stability is lost when 
the radius becomes in the order of the lattice constant. We 
consider the two-dimensional ferromagnet with up spins. The 
spin at the Skyrmion center is precisely oriented downward, 
which is topologically protected. The remnant of a Skyrmion 
burst is expected to be a single-spin flipped (SSF) point, which 
is a single down spin in the up-spin ferromagnet. It has a quan- 
tum diffusion dynamics governed by the Heisenberg Hamilto- 
nian. When the spin is S = h, solving this problem exactly, 
we find that the down spin hops about two-dimensional lattice 
points without canting. Namely, a SSF point can be regarded 
as a two-dimensional continuous-time quantum walker. We 
have verified manifestations of a quantum walk which differ- 
entiate from a classical random walk: The probability density 
takes the maxima at moving fronts and there appear oscilla- 
tions between moving fronts. We then analyze an interference 
of two SSF points after two Skyrmion bursts. The diffusion 
dynamics is exactly solvable also for S = 1. When the spin 
is higher (S > i), a SSF point breaks into 25 quantum walk- 
ers. This is a new type of quantum walk, which we may call a 
multiple quantum walk. We are able to create quantum walk- 
ers experimentally at any points in a magnetic thin film, first 



by creating Skyrmions sequentially and then by letting them 
collapse simultaneously. 

A quantum walk is a quantum analogue of a classical ran- 
dom The quantum walk corresponds to the tun- 
nelling of quantum particles into several possible sites, gener- 
ating large coherent superposition states and allowing massive 
parallelism in exploring multiple trajectories. The quantum 
walk is expected to have implications for various fields, for 
instance, as a primitive for universal quantum computing and 
systematic quantum algorithm engineering. Recently quan- 
tum walks have been experimentally demonstrated using nu- 
clear magnetic resonance[8], trapped ions|j9 10], photons in 
fibre optics 11 1:1 and waveguides ! 12fl . Our work presents an 
additional example of quantum walks in magnetic thin films. 

Giant Skyrmion: We use the classical spin field of unit 
length, n = (n x , n y ,n z ), to describe the spin texture whose 
scale is much larger than the lattice constant a = 0.3[nm]. 
The Hamiltonian consists of the anisotropic nonlinear sigma 
term Hj, the DDI term Ho and the Zeeman term Hz, whose 
continuous versions read 
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with the exchange energy T, the single-ion anisotropy con- 
stant £, the DDI strength f2, and the Zeeman energy Az. 
The ground state is the spin-polarized homogeneous state, 
n = (0, 0, 1), under the field perpendicular to the plane. 
The simplest cylindrical symmetric spin texture is 



cos( 



n x {x) = -yjl-a 2 (r) 
n v (x) = -y/l- <7 2 (Y) sin 0, n z (x) = a(r), (4) 

with 9 the azimuthal angle. The function er(r) should satisfy 
the boundary conditions a(r) = — 1 at r = to avoid the 
multivalueness, and er(r) — > 1 as r — > oo to approach the 
ground state. The spin texture (@J has the Pontryagin number 
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Qsky = 1, and hence it describes a Skyrmion. In particular, 
there exists a disk-like spin texture(magnetic bubble domain) 
with a domain wall at r = R. By changing variable as <j(r) 
'. an ii r (r), the equation of motion without the DDI term reads 
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We solve it by setting |tanhr| = 1 in the right-hand side, 



<r(r) = tanh 
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with Ii (x) the modified Bessel function. This is almost the 
exact solution outside the domain wall. Actually it gives also 
an excellent approximation to the domain wall. The excitation 
energy is given by [4] 

E sky (R) = ^^-Q[Rhx^-R] + 7T^A z + A7rT, (7) 
£ dp or 

where the ground-state energy is subtracted, and dp is the 
thickness of the film. The Skyrmion radius R is determined 
by minimizing E{R) with respect to R. Since R is as large as 
1/im for typical sample parameters, we have called it a giant 
Skyrmion||4j. The present Skyrmion is stabilized dynamically 
by the competition among the DDI, the Zeeman effect and the 
anisotropy of the film. The number of down spins is of the 
order (R/a) 2 , which is about 10 7 for R = 1/im. 

According to a method of collective coordinate, the dynam- 
ics of the Skyrmion radius is governed by the time-dependent 
Ginzburg-Landau equation, 
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dR = L dE sky (R) 
dt dR 

(8) 

where L is the Onsager's constant. The Skyrmion is stable 
when R is a constant. It becomes unstable when the magnetic 
field increases beyond a certain critical field. The DDI can 
stabilize the Skyrmion spin texture no longer. In this case, 
solving ([8]), we find the Skyrmion radius to change as 
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where Rq is the initial radius. It shrinks to zero (R 
a finite time interval, 
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It follows from © that, as R — > 0, u (r) = 1 for r > and 
er (r) = — 1 for r = 0. The spin texture (0 becomes 

n(0) = (0,0,1), n(x) = (1,0,0) fora^O. (11) 

The Skyrmion collapses into a singular point, which is a SSF 
point. 



Skyrmion burst: As the Skyrmion collapses, it emits spin 
waves: About 10 7 down spins are flipped upward, and the 
Skyrmion excitation energy E sky (R ) is carried away, leaving 
only the energy £"ssf of a SSF point. 

To describe spin waves we parametrize the spin field as 



n x (x) = a(x), n y (x) = \Jl - a 2 (x) siwd(x), 
n z {x) = \J\ - a 2 {x) cos#(x), 



(12) 



so that the ground state is reached by setting a(x) — d(x) = 
0. Note that this parametrization is different from that of the 
Skyrmion @. The Lagrangian density is given by 
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where S is the spin per atom. The Hamiltonian consists of 
the nonlinear sigma model Hz and the Zeeman term Hz- We 
may neglect the DDI Hjj to discuss spin waves. Substitut- 
ing the above configuration into the Hamiltonian density and 
taking the leading order terms in a and §, we obtain 



(v^r] + (i + %)(- 2 +n d4) 



The Euler-Lagrange equations are easily written down. There 
exists the outgoing propagating wave solution, 

cr(r,t) = Jq (kr) sinwi — Nq (kr) cos uit, (15a) 
§(r,t) = J (kr) cosujt + N (kr) sinujt, (15b) 

up to an appropriate initial condition, where the dispersion 
relation reads 



huj = -k 2 
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The asymptotic behavior is a (r, t) = y/2/nkr sin(fcr — wt), 
etc., for r ^> R. It carries away spins and energies from a 
Skyrmion. 

Quantum walk: The Skyrmion collapses into a singular 
point in the ferromagnet. The remnant is a SSF point in a fer- 
romagnet. The Skyrmion can be no longer a classical object 
but a quantum object. 

Our concern is about the dynamics of a SSF point after the 
skyrmion collapse. We show that the dynamics is described 
as a two-dimensional continuous-time quantum walk as gov- 
erned by the anisotropic Heisenberg model. 

The Heisenberg model is expressed as 



4^ 



(16) 

from which the nonlinear sigma model ([1} follows together 
withT = (l/2)S' 2 Jand^ 2 = AD/ J a 2 . The Zeeman term 
is H z = -(A z /S) J2i The DDI is irrelevant. Here, the 
index i runs over the two-dimensional lattice points, i G Z 2 . 
The ground state \g) is the up-spin polarized state defined by 



S?\g) = S\g), S+\g) = 0. 



(17) 
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FIG. 1: (Color online) (a) The probability dentity at t = 5/ J and 
t = 10/ J. The horizontal axes are m x and m y . 



FIG. 2: (Color online) (a) Time evolution of the probability dentity 
Px (t) • The horizontal axes are the time t and the site m x . 



A SSF point is generated at m = (m„ m y ) G Z 2 by applying 
(S 1 ") 25 to the ferromagnetic ground state, 



M = (S") 2S \g) 



(18) 



where only the spin at the lattice site m is pointed downward. 

First we study the case S — i in details, since the problem 
is exactly solvable. It is straightforward to show that 
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(19) 



The energy of a SSF point is -Essf = (m\H\m), which is 
much smaller than the energy E&y(Ro) of a Skyrmion. 

The dynamics of a SSF point is governed by the 
Schrodinger equation iHdijj/dt — Hip. The time evolution 
reads \ip (t)) — e~ l n H |0), where |0) denotes the initial state 
containing a SSF point at site m = (0, 0). By applying the 
Hamiltonian to the state containing one SSF point, the point 
remains at the same site or is shifted to one of the neighboring 
sites. The resulting state is a coherent superposition of them. 
It is notable that the spin is strictly parallel to the z-axis: It 
never cants. The down spin diffuses by time evolution in this 
way. We can regard a SSF point as a quantum walker. 

The quantum diffusion process of a SSF point is governed 
by a continuous-time quantum walk. Continuous-time quan- 
tum walk in one dimension is studied in the refs. lfl3l [Till . 
The process is a two-dimensional extension of continuous- 
time quantum walk. In ref . II 1 311 . the probability is obtained 
by the direct calculation of infinite multiplication of matrix. 
In ref. [14], Laplace transformation is used. We present a new 
derivation using the generating function, which is an easy way 
to generalize to a higher dimension. 

Expanding the exponential, we write the time evolution as 



e-^ H \0) = J^^t (m) \m). 



(20) 



Based on a combinatorial method known in the random walk 
theory, we can show that the coefficient £ t (m) is determined 



from the following generating function, 
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where m — (m x , m y ). It is determined in a closed form, 

C t = e it(i+^)i\m x \+\ mv \ j^Jt)^ ( Jt ) , (22) 
with the use of the generating function of the Bessel function, 

= £>./„(*). 



exp 



(23) 



The probability at site m is given by Pt(m) = |£ t (m)| 2 . The 
total probability is conserved, Y^m ^* ( m ) = as ^ s eas ily 
checked based on the formula J2 n ( z ) = 

We have thus solved the diffusion problem analytically. We 
show the probability density as a function of site m in FigQ] 
The probability density takes the maximum value not at the 
center but at the fronts. Inside the maximum values, the prob- 
ability density exhibits an oscillatory behavior in the scale of 
lattice constant. They are characteristic behavior of a quantum 
walk. This is highly contrasted to that of a classical random 
walk, where the probability density is Gaussian with the max- 
imum value taking at the center and not oscillating. 

The probability density is factorized as a direct prod- 
uct of the probabilities along the x-axis and the y-axis, 
P t (m x ,m y ) = P t (m x )P t (m y ), with 



Pt{m x ) = \J\ mx \ (Jt)\ 



(24) 



We illustrate the time evolution of P t (m x ) in Figf2] Long 
after the Skyrmion burst, it behaves asymptotically as 



P t {m x ) ~ —-cos 2 

TTJt 



(25) 



The velocity of the front propagation is given by v — 2J/ir. 
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FIG. 3: (Color online) (a) Time evolution of the probability dentity 
Px (t). The horizontal axes are the time t and the site m x . 

We study interference effects of two Skyrmions bursts. We 
generate two Skyrmions in a thin ferromagnetic film. By ap- 
plying large magnetic field, two Skyrmions collapse simulta- 
neously and two SSF points are generated. They are quantum 
mechanical objects, and the probability amplitude is given by 
a linear superposition of each SSF points. The time evolution 
of these two SSF points are shown in Figj3] An interference 
can be seen between two SSF points. This is a manifestation 
of quantum mechanical properties of a quantum walk. 

Multiple quantum walk: We discuss how the scheme is gen- 
eralized to the system where the spin is higher than 5 = |. 
In so doing, it is adequate to recapitulate the diffusion process 
for 5 = | from a slightly different view point. For 5 = \, 
a SSF point at site m is described by |m) = S^lg). When 
the Hamiltonian ( fl6b acts on this state, the operator 5~ re- 
mains at the same site or is shifted to one of the neighboring 
sites, as implied by (T% . Namely, the operator 5~ itself can 
be regarded as a quantum walker. 

When the spin is 5, by generalizing the above picture, there 
are 25 quantum walkers at the initial SSF point, and then they 
diffuse. The state is described by 

\m\m 2 ,--- ,m 2S ) = S ml S m2 ---S m2S \g) 1 (26) 

implying that a quantum walker is present at sites m\ The 
number of 5~ is independent of time due to the conservation 
of the total S z . Indeed, acting the Hamiltonian H to the state 
(|26] |. we find a coherent superposition of these states. The 
initial SSF point is given by d26b with m % — (0, 0) for all i. 

When 5 = 1, there are two quantum walkers. It is easy 
to calculate H\m , m 2 ) explicitly. The time evolution is also 
calculable. Since the Clebsch-Gordan coefficient is constant, 
we obtain a concise formula, 

€ t {m\m 2 ) =exp[it(j- A z - DS m i m 2)]£ t (m 1 )(r t (m 2 ). 

(27) 

The probability density is factorizable, Pt(m 1 ,m 2 ) = 
P t (m 1 )P i (m 2 ), with (E). 

However, since the Clebsch-Gordan coefficient is not con- 
stant for 5 > §, 

S- \q) = V 5/2(5/2 + 1) + 1) \q - 1) , (28) 



the analysis is not simple. Each walker does not diffuse in- 
dependently but interacts each other. Nevertheless, to get a 
rough picture, we dare to approximate (5~) 9 |(j) = C\q — 1) 
with a certain constant C. Then, we find 

2S 

P t (m\m 2 ,--- ,m 2S ) = [[^K). (29) 

q=l 

We may regard the diffusion process as a diffusion of 25 in- 
dependent walkers in this approximation. 

We have studied the diffusion process at zero temperature. 
One might wonder if it would survives spin- wave excitations. 
There exist no problem because spin-wave excitations have a 
large gap induced by the anisotropy and the external magnetic 
field. 

We have obtained an analytical solution of the time evo- 
lution dynamics of the SSF point generated by a Skyrmion 
burst. Its quantum diffusion process is a spreading oscillatory 
propagating wave and shows an interference. This is highly 
contrasted compared to a classical diffusion process, where 
the dynamics is described by a Gaussian. Our system provide 
a new approach to investigate properties of quantum walk. 
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